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SUMMARY 


An analytical analysis of conditions producing pilot induced oscillations 
(PIO's) was made for the Space Shuttle orbiter in a landing^approach configuration 
for the task of nulling the elevation angle of the line of sight to a target vehicle. 
The analysis yielded a value of PIO frequency and a value for the amount of total- 
system time delay (pilot + control system) that can be tolerated before instability 
results. Calculations were performed showing the effect of varying the range to the 
target and of varying the handling qualities of the orbiter vehicle. Comparisons 
were made of the analytical predictions and the simulation results obtained using the 
Langley Visual/Motion Simulator. Similar trends for PIO frequency were obtained for 
calculated and experimental values due to changes in vehicle handling qualities; how- 
ever, different trends were noted for changes in target range. Calculated results 
showed an increase in PIO frequency with a decrease in range, whereas the simulation 
results were invariant with range. In addition, calculated results showed a reduc- 
tion in the tolerable amount of time delay with a reduction in range which was sup- 
ported by the experimental results. 


INTRODUCTION 

A recent simulation study (ref. 1) was made using the Langley Visual/Motion 
Simulator to determine the effect of varying the amount of control-system time delay 
on the occurrence of pilot induced oscillations (PIO's) and on pilot tracking per- 
formance of a Space Shuttle orbiter in a landing-approach conf iguration. The pilot's 
task was to track a target airplane that performed a maneuver in altitude only. The 
target was initially offset laterally so that both orbiter pitch and roll hand- 
controller inputs were required. PIO's were encountered in both the longitudinal and 
lateral modes of motion when 250 milliseconds of additional delay were inserted in 
the simulation. 

The purpose of this paper is to present a PIO analytical analysis for the longi- 
tudinal short-period motion of the air-to-air tracking task that was simulated in 
reference 1. The analysis permits determining a value for the PIO frequency and a 
value of total-system time delay (pilot + control system) beyond which the system 
goes unstable. Several factors affecting the PIO are examined, such as target range 
and simulated vehicle handling qualities. Comparisons are made of the simulation 
results not included in reference 1 and the analytical calculations. 


SYMBOLS 

Numerical values are given for some quantities in both the International System 
of Units (SI) and in U.S. Customary Units for convenience. Measurements and calcula- 
tions were made in U.S. Customary Units. The effective stability derivatives used 
herein are referenced to a system of axes with the origin at the vehicle center of 
gra vi ty . 

A,B,C letter designates transfer function listed in table AI of appendix A 

a,b,c coefficients of polynomial listed in appendix B 
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vehicle mass, kg (slugs) 

period of oscillation, sec 

orbiter pitching angular rate, rad/sec 

2 2 

dynamic pressure, N/m (Ib/ft ) 

Laplace variable 

2 2 

wing area, m (ft ) 

orbiter velocity, m/sec (ft/sec) or knots 
orbiter weight, N (lb) 
pilot transfer function 

combined pilot and time-delay transfer function 
time-delay transfer function 
angle of attack, rad 
flight-path angle, rad 
determinant defined by equation (A3) 
hand-controller input after scaling, rad 
vertical tracking error (h^ - h) , m (ft) 
line-of-sight angle of pilot (fig- 1), rad 
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C damping ratio of longitudinal short-period mode 

9 pitch angle, rad or deg 

0 angle of target from orbiter measured positive above horizontal (see 

fig. 1 ) , rad 

X total time delay { / sec 

T control- system time delay, sec 

a 

T pilot effective time delay, sec 

e 

T maximum tolerable time delay before instability results, sec 

PIO ^ 

0 ) frequency of oscillation, rad/sec 

0 )^ undamped natural frequency of longitudinal short-period mode, rad/sec 

frequency of pilot induced oscillation, rad/sec 
Abbreviations: 

ALT approach and landing tests 

c.g. center of gravity 

PIO pilot induced oscillation 

VMS visual/motion simulator 

A dot over a symbol indicates the first derivative with respect to time. 


BACKGROUND 

Over the years longitudinal short-period PIO's have been encountered with a 
number of airplanes. A variety of different factors have been identified as the 
causative agents in generating these PIO's. Several theories have been advanced to 
try and predict PIO tendencies. The theory in reference 2 is based on pitch atti- 
tude tracking. Reference 3 provides another theory and also gives a bibliography of 
PIO papers. In reference 3 it is postulated that PIO's are generated by the pilot 
responding to normal acceleration cues. In the Shuttle orbiter tracking study on the 
Langley VMS it was believed that the visual task was the main driver in the simula- 
tion. Both vehicle pitching and altitude motions were involved. In addition, 
there were several nonlinear elements in the simulation. Since the theories of ref- 
erences 2 and 3 were not applicable to the task simulated, the present paper presents 
a PIO analytical analysis tailored specifically to the line-of- sight tracking task 
studied. The purpose was not to match prediction and experiment exactly, since the 
simulation involved nonlinear elements and both longitudinal and lateral modes of 
motion, but instead to indicate the trend linear theory would predict when varying 
the same parameters as varied in the simulation. 
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LINE-OF-SIGHT TRACKING TASK 


The six-degree-of-f reedom simulator study of reference 1 employed an air-to-air 
tracking task. For each run, the pilot was to visually track a target airplane as 
it varied slowly in altitude. Target altitude was driven by a low-frequency cosine 
wave. For each run only 1/2 cycle of target motion was employed so that the 
resulting piloting task would approximate a flare to the landing maneuver. For con- 
venience, the orbiter was in a level-flight trim condition at the start of each run. 
Also, the orbiter was at the same altitude but displaced laterally from the target. 

A reticle was supplied in the subject's visual field that was held fixed with respect 
to the orbiter body axes- The pilot's task was to keep the cross hairs on the target 
both in elevation and azimuth. 


The analytical analysis in this paper is concerned only with the orbiter motion 
in the plane of symmetry. For this situation the pilot's task was to null the eleva- 
tion line-of-sight angle sq. The geometry of the situation is depicted in figure 1. 
In addition to the target aircraft and reticle, the horizon was also displayed to the 
pilot. 


An examination of figure 1 shows the line of sight to be from the orbiter c.g. 
to the target c.g. This simple geometry is considered herein. The offset of the 
pilot's eye, forward of the orbiter c.g., and the location of the tail pipe, aft of 
the target c.g., complicate the geometry. Such arrangements, however, can be con- 
sidered. The horizon display is of use to the pilot in controlling bank angle, which 
is, of course, a lateral task. As indicated in figure 1, the pilot's longitudinal 
task involves only a single input, namely eg. This geometry yields the multiloop 
system shown in block diagram form in figure 2. 

The various elements of the block diagram are labeled and most are self- 
explanatory. The block marked time delay represents a transport time delay existing 
in the control system and was the primary variable in the simulation study of refer- 
ence 1. One of the basic assumptions made in analyzing the block diagram of figure 2 
is that the orbiter speed V remain constant. Another assumption is that the 
longitudinal distance between orbiter and target also remain constant; that is, 
is a constant. As a consequence, the dynamics block is reduced to only two degrees 
of freedom, which is the number required to represent the orbiter longitudinal, open- 
loop , short-period characteristics . 


The pilot's task, as depicted in figure 1 and as modeled in figure 2, is defined 
strictly as a line-of-sight tracking task. An alternate approach that is sometimes 
used to model this tracking task is to include an additional pilot-model block in the 


outer loop (see ref. 4) following the signal 


V 


in figure 2. 


Use of such an 


additional pilot block was omitted in the present study because of the descriptions 
given by the astronauts and pilots of the task performed in the simulator. As 
pointed out in appendix C of reference 1 , target aspect cues were seldom used by the 
subjects. All subjects indicated that keeping the cross hairs on the center of the 
tail pipe was the tracking goal in the simulation. Also, the subjects indicated that 
for the longitudinal task the horizon was of little use. 
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PIO CONDITION 


The PIO condition arises from the circulation of the pilot’s input signals in 
the various loops of the system. For a PIO to exist for the system of figure 2, the 
use of the low-frequency forcing function h^ is not necessary. Therefore, in a PIO 
analysis the input h^ can be set equal to zero. Also, in most PIO situations the 
pilot's input signals seem to have a wave form that closely approximates a sinusoid. 
This means that an analysis, if the system is linear, need only be concerned with the 
system frequency response. Finally, the PIO condition is defined to occur at that 
frequency corresponding to the largest value of time delay that can be tolerated 
before the system goes unstable. 


ANALYTICAL PROCEDURE 

First the characteristic equation for the system of figure 2 was determined. 

A mathematical model of the pilot was selected and inserted into this equation. By 
letting S = i(o, the stability boundaries were calculated in terms of two system 
parameters, pilot gain Kp and total-system time delay t. From a parameter-plane 
plot of Kp versus t, the largest value of t before unstable operation and the 
corresponding value of Kp were extracted. For the analysis herein, these values 
correspond to the PIO condition. Note that the value of co associated with these 
values of Kp and t is labeled the PIO frequency. 


ANALYSIS 


The equations of longitudinal motion for the vehicle dynamics block of figure 2 

were used to obtain the transfer functions and Details of the development 

o o 

e e 

are given in appendix A. 
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given, two transfer functions , one representing control- 


e e 

system time delay and the other representing the human pilot, need to be defined. 
Recall that the transfer function which represents a pure time delay can be written 
in Laplace notation as 


Y 


TD 



S 


( 1 ) 


For the pilot model in the inner loop of figure 2, the transfer function selected is 


Y 


P 


K 



( 2 ) 


Several choices for a pilot model exist in the literature, from the "synchronous 
pilot" represented by a simple gain Kp (see ref. 5 and background review in 
ref. 3) , to a very sophisticated model (see ref. 6) that contains lead and lag equal- 
ization terms, a time delay, and second-order systems for neuromuscular effects. The 
pilot model represented by equation (2) discards the lead and lag equalization terms 
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of the sophisticated model from consideration for the PIO condition; however, it 
retains the time delay and the lag effects of the neuromuscular system through the 
use of an effective time delay 

Since time delay appears in both equations (1) and (2), it was convenient math- 
ematically to combine these two transfer functions as depicted in figure 3(a). Thus 


^PTD = 


-ts 


(3) 


where 


T = T + X 
e a 


(4) 


Once the various transfer functions of figure 2 are specified, an analysis of 
the closed loop can be undertaken and the system characteristic equation can be 
derived. Details of the derivation are presented in appendix B. The characteristic 
equation thus obtained can be written as 


A + Y 


PTD 



(5) 


where A, B, and A are defined in table AI. Substituting equation (3) into equa- 
tion (5) yields 


A + Kpe 


-xSj 


AS + 


VB 


= 0 


(6) 


This equation involves the two unknown parameters Kp and x. To obtain stability 
boundaries, substitute S = ico into equation (6) and determine the values for Kp 
and X which satisfy the resulting equation. 

In solving an expression such as equation (6) that involves two unknown param- 
eters, the general technique is to rewrite the single equation as two equations. One 
of the new equations deals with the real part of the original expression and the 
other deals with the imaginary part. These two equations with two unknowns can then 
be solved simultaneously to determine values of the unknown parameters. However, in 
the particular case of equation (6), a different approach is necessary since it 
involves an exponential term. This different approach solves for the amplitude and 
phase angle of Yp^p* Evaluation of Kp is then possible since it is associated 
only with the amplitude, and similarly, evaluation of x is possible since it is 
associated only with the phase angle. (See ref. 7.) 


An illustrative parameter-plane plot of Kp versus x is given in the sketch 
in figure 3(b). The curve is a stability boundary that divides the parameter plane 
into stable and unstable regions. In the sketch the maximum value of x is labeled 
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as the PIO condition. This is the largest value of time delay that can be tolerated 
before the system goes unstable. With this value of t there exists a single 
corresponding value of Kp. Associated with these two values is a frequency called 
the PIO frequency 

In figure 3(c) additional detail is presented that helps illustrate the reason- 
ing behind the selection of the PIO condition. Recall that t is the sum of two 
delays, and Since the pilot's effective time delay is unchanging, and 

therefore a constant, there must exist an inaccessible region in the parameter plane 
as shown in the sketch. That is, x = Tq is the minimum possible delay value. Of 
course, if x - Xq the amount of control-system time delay x^ is zero. The 
range of Kp represented by a vertical line at this value of delay and contained 
within the stability boundary is quite large. From an examination of figure 3(c), it 
is apparent that adding an increment of control-system time delay x^ results in a 
decrease in the range of possible pilot gain for stable system operation. Continuing 
to increase x^ constrains the stable range of pilot gain even more. Finally, there 
exists a single value of delay Xp;£0 beyond which the system is unstable no matter 
what the value of pilot gain. Thus, it is apparent that for any value of x less 
than xpiQ, the pilot has the ability if he finds that he is operating in an unstable 
region of the parameter plane to achieve stable operation by simply altering his 
gain. For the PIO condition, the pilot can no longer affect system stability by 
adjusting his gain. 

It is, of course, recognized that any point along the boundary defines a stabil- 
ity situation that can be labeled as a PIO condition. The use of rationale other 
than that just described could lead to selecting any number of different points along 
the boundary curve as PIO conditions. As a consequence, different values of pilot 
gain, time delay, and frequency would result. For this paper, constraining the 
pilot's control action in the manner described previously seems an appropriate choice 
for the tracking task simulated. 

It is of interest to note that in the simulation study of reference 1, the hand- 
controller characteristics were nonlinear. For this analysis, the hand-controller 
characteristics are assumed linear. Controller characteristics enter the pilot 
transfer function as a single constant gain and thus are contained in Kp. 


CASES CONSIDERED 

The cases calculated were chosen to provide a range of values that would 
encompass the test conditions of a set of simulator runs obtained as a sequel to the 
study of reference 1 . The analytical calculations used the stability derivatives and 
physical characteristics given in table I. These values are identical to those used 
in reference 1 for the Shuttle orbiter. Table II provides expressions for converting 
the nondimens ional derivatives of table I to the dimensional form for use in the 
mathematical development herein. 

The cases calculated are concerned with examining the effect of range and the 
effect of changing vehicle handling qualities. Calculations were performed for six 
equally spaced values of range, from 30.48 m (100 ft) to 182.88 m (600 ft), for each 
different set of parameter values considered. Vehicle handling qualities were varied 
by using the four different sets of parameter values given in table III. Three 
sets of dimensional parameters were obtained by simply increasing or decreasing the 
orbiter speed from the value used in reference 1 (basic case) . The remaining set was 
obtained by changing the single stability derivative C 2 ^ of the basic case. 
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RESULTS AND DISCUSSION 


General Remarks 

‘The analytical development presented herein provides the following two numerical 
results: 

1 . A value for the PIO frequency 

2. A value for control-system time delay beyond which the closed-loop system 

goes unstable 

The calculated results for these two entities are discussed in the next section 
entitled "Analytical Results." The effects of range and changing orbiter handling 
qualities are covered briefly. The purpose of the calculations was to provide values 
that correspond to simulator supplemental test data. Explanations and the reason for 
performing the supplemental simulator study are presented in the "Experimental 
Results" section. A comparison section discusses the trends obtained from theory and 
from experiment. 


Analytical Results 

Calculations were performed to determine the values of pilot gain Kp and time 
delay t as a function of frequency for six different values of range between 
orbiter and target for each of the four different handling-quality cases given in 
table III. Values of range from 30.48 m (100 ft) to 182.88 m (600 ft) were 

used. Parameter-plane plots of Kp versus x were constructed at each value of 
range and the PIO conditions were determined. Most of the parameter-plane plots were 
as sketched in figures 3 and 4(a). For the shorter ranges, however, a loop appeared 
as shown in figure 4(b) that enlarged in size as range was reduced. For these situa- 
tions the PIO condition occurred at the loop closure point. The PIO conditions 
obtained from all of the parameter-plane plots are given in figures 5 and 6 showing 
PIO frequency (x)pio time delay Xpjo ^ function of range between orbiter and 

target. 

All of the curves shown in figures 5 and 6 indicate that 0)pio increased with 
a reduction in range. In addition, as range was reduced the amount of time delay 
Xpio that could be tolerated in the system before instability resulted was also 
reduced. As a consequence, these results point out the dilemma facing a pilot per- 
forming the task of closing on a target. It is readily apparent that as range is 
reduced, some value of range will be encountered at which a PIO will occur and below 
which the pilot-vehicle system will go unstable. 

Changing vehicle handling qualities is accomplished herein in two different 

ways. One way is by changing the orbiter speed and the other way is by changing a 

single derivative Cy . Both ways were used in the supplemental simulator tests. 

a 

The values of V and Cy used here are the same as those used in the simulator 

runs. Results for the effect of handling qualities on PIO conditions of frequency 
0 )pjQ and time delay xpxo assessed using the different curves in figures 5 

and 6. Increasing either speed or the magnitude of ^2^ improves the vehicle 
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handling qualities. At any given value of range, an increase in V or the magnitude 
of increases ^^10 reduces xpio* 

It is worth noting that the total amount of time delay in the system TpjQ is 
used as the ordinate in figures 5 and 6. From equation (4) recall that the total 
time- del ay value consists of the sum of two parts, one part due to the pilot and 

the other part due to the control system To ascertain the amount of time delay 

at the PIO condition due to the control system requires that an estimate be made of 
the pilot delay term x^ and the difference determined. A typical value of x^ of 
about 200 msec is believed to correspond to the normal experience of a pilot perform- 
ing a conventional tracking task. A value of Xq of 100 msec can be considered as a 
lower limit or bound of pilot effective time delay when the pilot is represented by 
the simple transfer function of equation (2) and when he is pushed to his operating 
limit. Such a case might occur during a PIO. Simply reducing the ordinate values of 
Xpio in figures 5 and 6 by a value of Xg of 100 msec provides the largest values 
of control-system time delay x^^ that can be tolerated before the pilot-vehicle 
system becomes unstable. 


Experimental Results 

The simulator study of reference 1 using two astronauts and two research pilots 
as test subjects gave periods for the longitudinal PI0*s of about 5 sec for all test 
subjects. The period of the PIO obtained on landing number 5 of the Shuttle orbiter 
at Dryden Flight Center was about 2.5 sec. Because of this difference, a few 
additional tests were made in the Langley Visual/Motion Simulator using one of the 
research pilots in the original study to see if a 2. 5- sec period could be obtained 
with the simulator. 

In the original simulation, the PI0*s were found to occur near the end of the 
run. For those tests the range between the target vehicle and the orbiter was varied 
linearly with time. Initial range was 182.88 m (600 ft) and final range was 91.44 m 
(300 ft). Run time was 60 sec. At run termination the target wing span in the 
visual scene matched the width of the horizontal bar of the reticle. This scheme was 
used to induce the pilot to increase his gains as he normally does during the landing 
flare. Because the PIO’s occurred at the closer ranges, some supplemental runs hold- 
ing range constant were made at even closer ranges. Runs were made at 91.44 m 
(300 ft), 60.96 m (200 ft), and 30.48 m (100 ft) to determine the effect of range on 
PIO period. The runs were made with various amounts of time delay present. PIO’s 
did not occur for all of the runs made; however, sufficient oscillations occurred in 
the a and 0 time-history traces to determine a value for the period due to con- 
trol inputs. The results of these tests indicate that the period of the oscillatory 
motions was invariant with range. Some additional runs were then made for different 
vehicle speeds. The original simulation setup was for a 300-knot (500 ft/sec) 

Shuttle orbiter. Two additional speeds were tested with one greater and one less 
than the original 300-knot speed. At each speed several runs were made at the dif- 
ferent ranges. Figure 7 presents the results of these tests. Note that changing 
speed changes the handling qualities of the vehicle as shown in figure 8 (ref. 8). 

The results show that for the higher speed configurations (better handling quali- 
ties), the period of the oscillatory motions was reduced to about 3 sec. The period, 
however, remained invariant with range for the different speeds. 

Increasing speed affects both the longitudinal and lateral handling character- 
istics. Earlier tests have shown that an improvement in one channel usually results 
in an improvement in the other channel. To eliminate the possible influence of an 
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was changed froin -3 ,45 to -6.90 for 


easier lateral task, one single derivative 

the same 300-knot speed condition to give only a change in vehicle longitudinal char- 
acteristics. Several runs were made at the same three ranges. The results are shown 
in figure 9. The data for the Shuttle vehicle from figure 7 are shown for compari- 
son. Periods of the oscillatory motions in a and 0 were 3 sec for the modified 
configuration. Figure 10 (ref. 8) locates this modified configuration on the han- 
dling qualities chart. 

Figures 11 and 12 show the time-history records of several of the runs made. 

Note that in figure 11 the pilot tried to keep his gains high and, consequently, at 
about 45 sec the motion base drive signals exceeded some preset limitations and the 
motion base automatically reset to the neutral position. In figure 12 the pilot 
simply stopped making inputs at about 45 sec. In both figures, eg is the line-of- 
sight angle and 6^ is the hand-controller input behind the quadratic shaper that 
goes directly into the equations of motion. 

The basic conclusions from these supplemental tests were that (1) the PIO fre- 
quency is invariant with range, (2) PIO frequency increases with an increase in speed 
(improved handling qualities), and (3) PIO frequency increases with an increase in 
C 2 ^ (improved handling qualities). It is important to note that the main purpose of 

these supplemental tests was to establish the frequency of the PIO. No effort was 
made to assess the amount of control-system time delay necessary for PIO to occur at 
any given test condition of range and speed. 

It should be pointed out that the target vertical movement during the run and 
the orbiter lateral offset at run initiation were arbitrarily reduced as the range 
was reduced. These changes were introduced in order that the elevation and 
azimuth visual angles presented to the subject during the run remained at about the 
same order of magnitude. Also, for the data presented herein, the offset of the 
pilot’s eye from the orbiter c.g. was eliminated. Finally, in performing the sup- 
plemental tests the same amount of time delay was inserted in both the pitch and roll 
control channels immediately following the pilot’s control input signals. The set 
of time-delay values used were 0, 4, 8, and 12 units of time delay as used in refer- 
ence 1 where each unit represents the update interval of the digital computer used. 
(Each unit represents 31.25 msec.) Each of these values has been increased by 
50 msec to account for the delay present in the digital computer from input to output 
and for the delay present in the visual scene generation equipment. (See ref. 1.) 
These are the values listed as in figures 7 to 12. 


Comparison of Results 

A comparison of experimental results and analytical predictions for the PIO con- 
dition is given in figure 13. Separate comparisons are presented for each of the 
four different handling-quality cases tested. An examination of the four frequency 
plots shows, as pointed out earlier in the text, the invariance with range of the 
experimental data and the increase in frequency with decreasing range of the analyt- 
ical results. Of additional interest is the fact that the predictions give a larger 
value of PIO frequency than was obtained experimentally for all test conditions. 

Comparisons of the PIO condition for time delay are also given in figure 13. 

The total value of time delay Tpjo is presented. In order to make the compari- 
sons, a value of Xq had to be added to the experimental delay values used for the 


11 



simulator tests. A value for of 100 msec was used. Also, the simulator had 

about a SO^msec delay due to the digital computer and visual display employed. (See 
previous section and also ref. 1.) Therefore, the minimum value of x that could be 
tested was about 150 msec. Figure 13 shows several tests performed at this condi- 
tion. In examining the figure it is interesting to note that at the closest range 
tested even this minimum value of time delay exceeded the values of the calculated 
curve. In fact, for most tests in the simulator the values of time delay used 
exceeded those predicated as the PIO condition by the analytical analysis. Of course 
at the time the tests were made the analysis had not been initiated. Notice that 
shading of the time-delay symbols is used to distinguish those cases in which PIO*s 
occurred. When the experimental test points lie below and to the right of the bound- 
ary no PIO*s occurred. In the region of the boundary when the PIO’s occurred they 
occurred near the end of the simulation run. When the experimental points were far 
to the left or above the curve a definite instability occurred early in the run and 
the time-history traces of the PIO were very erratic throughout. For such circum- 
stances the analytical analysis indicated that an unstable situation existed. What 
occurred in reality was that the pilot altered his tracking task and no longer 
tracked the target as postulated in the mathematical analysis. 

It is worth pointing out that the research pilot used for the supplemental tests 
and all of the test subjects of reference 1 indicated that the task in the simulator 
was a line-of-sight tracking task as modeled herein. It is possible, however, 
by altering the model to obtain theoretical values for copjo that are invariant with 
range. This can be accomplished by including an additional pilot gain constant 

in the outer loop of figure 2 with e — as the input. It is readily apparent 

that adjustment of the value of this additional pilot gain to offset the effects of 
changes in V ^ould provide results invariant with range. Note that this is 
equivalent to subjects changing their aim point to other locations inside the tail- 
pipe area rather than the center. Acceptance of an area as the aim point is believed 
to be what the pilot used as a practical goal in performing the task in the simulator 
at the closer ranges. The subject comments in appendix C of reference 1 support this 
altering of the tracking goal at the closer ranges. 

To facilitate examining the effect of vehicle handling qualities on frequency, 
figure 14 presents a comparison of analytical predictions and experimental results 
for a single value of range of 91.44 m (300 ft). The differences in magnitude 

between experimental and calculated values could be expected from the comparisons in 
figure 13. The trends of the results, however, with improving handling qualities 
either by increasing the vehicle's speed or by changing the single derivative 
were comparable. 

It is important to note that there are other possible reasons for differences 
existing between the simulation results and the analytical predictions, such as 

1 . The simulation involved six degrees of freedom and the tracking task required 
the pilot to make both pitch and roll control inputs for task accomplishment. For 
the supplemental tests the same amount of time delay was inserted following the pitch 
and roll control input signals from the hand controller. When the pilot experienced 
control difficulties longitudinally he nearly always had difficulties laterally. 

Thus, the experiment involves a multi-axis control task whereas the theory is for a 
single-axis task. 
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2. Motion cues were present in the simulator and not accounted for in the 
theory. 

3. A quadratic shaping circuit modifying the hand-controller output was present 
in the simulation and not in the theory. 

4. Subject awareness of potential resetting of the motion base to the neutral 
position because of overly aggressive control action could have resulted in an exper- 
imental task at the closer ranges different from the theory. 

There is one other comparison of PIO frequency that needs to be considered, and 
that is a comparison between the PIO periods obtained in the simulation of the basic 

case (P = 5 sec) and in the full-scale orbiter flight tests at Edwards Air Force Base 

(P « 2.5 sec). One possible reason for the difference is that the orbiter aerody- 
namics at the time of the PIO were not identical. The PIO of free flight number 5 
occurred just prior to touchdown. For this condition the vehicle was in ground 
effect. For the simulation, the orbiter was modeled at an altitude of about 457 m 
(1500 ft). Thus, some differences exist in the orbiter aerodynamic characteristics 
for these two cases which could be expected to affect the PIO periods. In addition, 

the piloting tasks were not the same. For the flight tests, the pilot was also con- 

cerned with modulating the flight velocity, executing a landing flare, arresting the 
sink rate, and achieving the desired attitude for touchdown. These additional con- 
cerns were missing in the simulation; nevertheless, it should be noted that the 
pilot’s task in free flight number 5 was to land at a particular spot on the runway. 
Thus, as the orbiter proceeded along the glide slope the pilot was performing a 
tracking task to achieve the touchdown point that is similar to the task modeled ana- 
lytically herein. The fact that a PIO occurred for the Shuttle orbiter on free 
flight number 5 should not be too surprising since the analysis herein shows that for 
a line-of-sight tracking task the amount of time delay that can be tolerated in the 
system for stable operation decreases as range to the target decreases. 


CONCLUDING REMARKS 

An analytical analysis of conditions producing pilot induced oscillations 
(PIO’s) was made for a Shuttle orbiter configuration for the task of nulling the 
elevation angle of the line of sight to a target vehicle. The analysis yielded a 
value of PIO frequency and a value for the amount of total-system time delay that 
can be tolerated before instability results. Calculations were performed showing 
the effect of varying the range to the target and of varying the handling qualities 
of the orbiter vehicle. Comparisons of the analytical predictions were made with 
results from a six-degree-of-f reedom simulation using the Langley Visual/Motion 
Simulator. 

Results of the analytical analysis indicated an increase in PIO frequency with a 
decrease in range, whereas the simulation results were invariant with range for the 
ranges tested. Predicted values of total-system time delay that could be tolerated 
before PIO occurs decreased with a decrease in range. For both calculated and 
experimental results, improving the vehicle handling qualities either by increasing 
the vehicle speed or by increasing the magnitude of the single stability derivative 
C^^ increased the frequency and reduced the tolerable value of time delay. 

The analytical analysis presented herein was made for the purpose of indicating 
what linear theory would predict for a visual line-of-sight tracking task. It was 
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felt that the results supplement the simulator data and perhaps offer some explana- 
tion for the variations obtained. It is recognized that a number of factors other 
than time delay can result in PIO*S/ such as control-system nonlinearities, rate 
limiting of control surfaces, acceleration effects, and others. This study also 
indicates that PIO can be task dependent. Thus, the results are not believed to 
be directly applicable to analyzing the longitudinal PIO obtained in free flight 
number 5 of the approach and landing tests (ALT) since the piloting tasks are not 
the same. It is felt, however, that a similar trend toward instability with decreas- 
ing range for a line-of-sight tracking task would apply. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 19, 1982 
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APPENDIX A 


EQUATIONS OF MOTION, SHORT-PERIOD CHARACTERISTICS, 

AND OPEN-LOOP TRANSFER FUNCTIONS 

The standard linearized longitudinal perturbation equations of motion are used 
for this mathematical development. They are further simplified by the following 
assumptions: 

1 . Constant airspeed 

2. Initial condition is straight and level flight 

3. Stability derivative M* is neglected 

Under these conditions the equations of motion, written in Laplace notation are 


(S + L ) a + (L - 1)q = -L. 5 
a q ^ 6 e 

e 


(A1) 


( -M ) a + ( S 
a 


M 

q 


)q = Mg 


6 

e 


e 


(A2) 


The longitudinal 
equations (A1) and (A2 


short-period mode characteristics can be easily obtained from 
) • For convenience 


A = 


(S + L ) 

OL 


("M ) 
a 


(L - 1) 

q 

(s - M ) 
q 


(A3) 


The characteristic equation can be obtained by expanding the determinant in equa- 
tion (A3) and setting the result equal to zero. Thus 

+ (L - M )S + M (L - 1) - L M =0 (A4) 

a q Ct q a q 

From this 


0)^ = M(L -D-LM 
n a q a q 


(A5) 


2Cw = L - M 
n a q 


(A6) 
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L - M - 

« q 

2fM (L - 1) - L 


CT 0 

^ Using equations (A1 ) and (A2) the open-loop transfer functions and 

-g— can be obtained directly as e e 


3 ^ = 
6 
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(s + L ) 
a 

(-M ) 
a 

{■\) 

i\) 

(s + L ) 

(L - 1) 

a 

q 

(-M ) 

(S - M ) 
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i-\) 
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- 1) 
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i\) 

(S - 

M ) 

q 

6 

e 

(S + L ) 
a 

(L 

q 

- 1) 


(-M ) 

(s - 

M ) 


Table AI lists these transfer functions in expanded form* 

The altitude transfer function can be obtained by recalling (see fig. 2) 

Y = 9 - a 


V 


Combining equations (A11 ) and (Al2) yields 


^ vl ■ 
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Rearranging equation (A13) gives 


6 



(A14) 


and finally 


h 

6 


V 

S \6 



(A15) 


With the appropriate substitution, equations (A14) and (A15) are expanded and given 
in table AI. 
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APPENDIX B 


CLOSED-LOOP ANALYSIS 

The following e;q)ressions were developed using the block diagrams of figures 2 
and 3(a). From the diagrams 


e 

V 



h 


Rewriting h as 



(B1) 


(B2) 


and inserting into equation (B1) yields 



(B3) 


Rearranging gives 




(B4) 


The characteristic equation for the system of figure 2 can be obtained by setting the 
denominator of equation (B4) equal to zero. Thus 



(B5) 


Rewrite equation (B5) as 


1 



(B6) 
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Note that 



e 

V 


1 



(B7) 


Substituting this identity into equation (B7) and rewriting gives 


1 + 


^0 . 
- — = 0 

0 t 


(B8) 


but 


^0 ^ 1 


(B9) 


where 


Y 


e 



(B10) 


Substituting equation (B9) into equation (B8) and rearranging gives 


1 + + 



0 


(B11) 


Now using appendix A and figure 3(a) rewrite equation (B10) as 


Y = — Y 

0 S A PTD 

Similarly 



(B12) 


(B13) 
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and 


^ = -5^ y 
^9 A 


(B14) 


Inserting equations (B12) and (B14) into equation (B11) yields 


1 + — ^ Y + Y =0 

S A PTD i 2 . PTD 
t S A 


(B15) 


Combining terms yields 


A + Y (as + — ^ = 0 

PTD I 


(B16) 


Substituting for A and B from table AI and expanding yields 


S^ A + Y (aS^ + bS + c) =0 
PTD ' ' 


(B17) 


The coefficients a, b, and c are given in the following table and the expression 
for A is given in table Alz 


Symbol 


Expression 


V 


X 5 0 

t e e 


M. L - ML- + — 
o a a 0 
e e t 


M. - L. M + (L 
6 o a o a 

e e e ^ 


- 1 ) 


V 


M. L - ml. 
\\ o a a o 
t \ e e 


Equation (B17) was used to compute the parameter-plane curves presented herein after 
first rearranging the equation to solve for and then substituting Kpe”“^^ 

for Yp^j^ as indicated in the text. 
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TABLE II.- DIMENSIONAL-DERIVATIVE CONVERSION 


FORMULA AND EQUATIONS FOR VEHICLE OPEN-LOOP 
LONGITUDINAL SHORT-PERIOD CHARACTERISTICS 


L 

a 


-qS 

C 

mV 2 

a 


L 

q 



mV 
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2V 


L 


6 

e 



mV 


5e 
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qS c 
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I 
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c 

m 
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qS c 

w _ c 

M = C 

q I m 2V 

y q 


M 


6 
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qS c 
w 



2 

CO = -(M + L M -ML) 

n a a q a q 


2Cco = L - M 
n ct q 
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TABLE III.- VALUES OF DIMENSIONAL DERIVATIVES AND VEHICLE OPEN-LOOP 


LONGITUDINAL SHORT-PERIOD CHARACTERISTICS 


Parameter 

Low speed 

Basic case 

High speed 

Modified 
basic case^ 

V, m/sec 

86.8680 

152.4000 

217.9320 

152.4000 

V, ft/sec 

285 

500 

715 

500 

V, knots (approx.) 

170 

300 

430 

300 

q, N/m^ 

4624.1077 

14 232.4067 

29 103.8472 

14 232.4067 

q, Ib/ft^ 

96.5765 

297.2500 

607.8465 

297.2500 

L 

a 

.5508 

.9664 

1.3819 

1.8727 

L 

q 

.1940 

.1940 

.1940 

.1940 


-.0917 

-.1609 

-.2301 

-.1609 

M 

a 

-.0399 

-.1229 

-.2512 

-.1229 

M 

q 

-1.8176 

-3.1887 

-4.5598 

-3. 1887 


.4665 

1.4359 

2.9363 

1.4359 

0)2 

n 

1.0333 

3.1805 

6.5037 

6.2619 

0) 

n 

1.0165 

1.7834 

2.5502 

2.5024 

C 

1.1649 

1.1649 

1.1649 

1.0113 

2 

n 

2.3682 

4.1551 

5.9417 

5.0614 


^Changed only to -6.90 from -3.45 of basic case. 
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Reference altitude 

Figure 1.- Tracking task geometry. (Positive values of variables are shown.) 
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(a) Section of closed-loop system. 



(b) PIO condition in parameter plane. 



(c) Arbitrary location in parameter plane. 
Figure 3.- Closed-loop analysis sketches. 
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170 


Hard to fly with large delays. 
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Expect lower values for 
smaller delays. 


(B) Data for small oscillation 
amplitudes only. 

For large amplitudes, 
periods of 4 to 4.5 sec occur. 

Time delay, 

Tg. msec 
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O 300 
A 425 


300 
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oversensitive 


Initial response fast, 
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V \3.5 / > 


Speed, 

knots 


■“Slow initially^ 
then oscillatory; 
tendency to 
overcontrol 


Very slow 
response 
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Damping ratio, ? 

Figure 8.- Effect of changing orbiter speed on handling qualities as shown 
on the longitudinal short-period handling qualities chart from reference 8 
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1 Time history with V = 300 knots, 
60.96 m (200 ft), and = 300 msec. 
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Time, sec 

Figure 12.- Time history with V = 430 knots, 
= 60.96 m (200 ft), and = 50 msec- 
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